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Modelling the Mechanics and Hydrodynamics of 
Swimming E. coli 

Jinglei Hu, a ’ b Mingcheng Yang, c Gerhard Gompper, a and Roland G. Winkler^ 


The swimming properties of an E. coli -type model bacterium are investigated by mesoscale hy¬ 
drodynamic simulations, combining molecular dynamics simulations of the bacterium with the 
multiparticle particle collision dynamics method for the embedding fluid. The bacterium is com¬ 
posed of a spherocylindrical body with attached helical flagella, built up from discrete particles 
for an efficient coupling with the fluid. We measure the hydrodynamic friction coefficients of the 
bacterium and find quantitative agreement with experimental results of swimming E. coli. The 
flow field of the bacterium shows a force-dipole-like pattern in the swimming plane and two vor¬ 
tices perpendicular to its swimming direction arising from counterrotation of the cell body and the 
flagella. By comparison with the flow field of a force dipole and rotlet dipole, we extract the force- 
dipole and rotlet-dipole strengths for the bacterium and find that counterrotation of the cell body 
and the flagella is essential for describing the near-field hydrodynamics of the bacterium. 


1 Introduction 

The bacterium E. coli is an example of a widely-studied class of 
motile microorganisms that exploit multiple helical flagella for 
locomotion. 1 Each flagellum is propelled by a reversible rotary 
motor anchored in the cell’s membrane. ® When all the motors 
rotate in the same direction, the flagella form a bundle and the 
bacterium swims forward; the so-called ‘run’ phase. It is inter¬ 
rupted by short periods of ‘tumble’ events, where the reversal of 
the motor-rotation direction of some flagella causes the associ¬ 
ated flagella to leave the bundle, thereby inducing erratic rotation 
of the cell body. 1 30 When the reversed motors switch back to 
their initial rotation direction, the bundle is reformed and the bac¬ 
terium swims in a new direction. The alternating runs and tum¬ 
bles allow the bacterium to efficiently execute a biased random 
walk toward favorable environments such as food-concentrated 
regions by adjusting run and tumble durations to the environ¬ 
mental conditions. 

Numerous experimental, theoretical, and simulation studies 
have been performed to unravel the physical aspects of bacteria 
locomotion.® These comprise aspects of the bacteria flagella such 
as their polymorphic transformationsEEGSHS] and bundle forma¬ 
tion. 15 E51 Moreover, bacterial propulsion properties have been 
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sir run-and-tumble dynamics J 16 l 24 l as well 
as the influence of hydrodynamic interactions on their motion ad¬ 
jacent to surface .E3E 35 In addition, the effects of external flows 
on the dynamical behaviors of bacteria suspension have been ad- 
dressed®MMI along with their rheological properties.E2B 1 


The complexity of the bundling and swimming processes, espe¬ 
cially near-field hydrodynamics, poses substantial challenges for 
an analytical description of bacteria locomotion. Here, mesoscale 
hydrodynamics simulations are particularly valuable to gain in¬ 
sight into the microscopic aspects of swimming, because they are 
able to bridge the large length- and time-scale differences be¬ 
tween the bacterium and fluid degrees of freedom J 42 l 43 i The mul¬ 
tiparticle collision dynamics (MPC) methodE3E3has proven to be 
very valuable for the studies of active systems JSEEmHI Specifi¬ 
cally, MPC has successfully been applied to elucidate synchroniza¬ 
tion between the flagella beating of nearby swimming sperm, EH 
as well as bundling of helical flagella of bacteria.EM 9 


Valuable theoretical insight into the swimming behavior of bac¬ 
teria is achieved by very simplified models. An ellipsoidal or 
spherical body is often combined with either a single effective 
flagellum attached at one pole, 13 l 35 l 55 l or with several flagella 
attached in a more or less random manner. ^ 24 ! 34 ! Thereby, hy¬ 
drodynamic interactions are typically taken into account by a hy¬ 
drodynamic tensor, 35 55 e g ^ Oseen or Rotne-Prager ten- 
sor.ES Such an approach is useful for a moderate number of hy- 
drodynamically interacting units, typically a single or very few 
flagella. In contrast, studies of large systems and collective phe¬ 
nomena involving multi-flagellated cells require a different ap¬ 
proach. We suggest to exploit the MPC method to account for 
fluid hydrodynamics. Moreover, a model for E. coli -type bacteria 
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that can reproduce experimental results and allows for quantita¬ 
tive predictions of their swimming properties is desirable; here, 
only a few modelling studies are available so far J— ^l 24 l 3 5[ 

In this article, we present a bacterium model which closely re¬ 
sembles the geometry, flagellar elastic properties, and rotary mo¬ 
tor torque of E. coli. By MPC simulations, we show that this model 
quantitatively reproduces the experimentally measured proper¬ 
ties of E. coli for both the hydrodynamic friction and the relation 
between the bacterial swimming speed and flagellar rotation rate. 
We find that the flow field near the model bacterium is rather 
complex with two spiral vortices arising from counterrotation of 
the cell body and flagella. At larger distances, the flow field dis¬ 
plays a dipole pattern comparable to that of swimming E. coli. We 
perform a detailed analysis of the flow by comparison with the 
force-dipole and rotlet-dipole approximations, which enables us 
to access the force-dipole and rotlet-dipole strengths, which are 
essential quantities for the description of bacterial swimming. 

The rest of the manuscript is structured as follows. Section [2] 
outlines the bacterium model and the MPC method for the fluid. 
Simulation results are presented in Sec. [3] and Sec. [^summaries 
our findings. Further technical details for the flagellum modelling 
are presented in App. [a] and the analysis of the flow field for 
periodic systems in App.jB] respectively. 

2 Simulation model and method 

2.1 Bacterium model 

The bacterium consists of a body and flagellar filaments, as shown 
in Fig. [l] which are composed of point-like particles of mass 
M = 10m. The cell body is represented by a spherocylinder of 
diameter d = 9a and length / b = 25 a, composed of 51 circular sec¬ 
tions of particles with a spacing of 0.5a, see Fig. |2fa). Here, m 
and a are mass and length units related to the MPC fluid, as de¬ 
scribed in Sec. |2.2| Both pole sections consist of a single particle 
only. In each of the other 49 sections, 60 particles are uniformly 
distributed along circles on the spherocylinder surface. To main¬ 
tain the shape of the body, nearest- and next-nearest-neighboring 
pairs of particles are bonded by a harmonic potential 

U bd =^K bd (r-r e )\ (1) 

where r and r e are the distance between the pair and the pre¬ 
ferred value, respectively. Moreover, two particles separated by 
10, 20, and 30 particles along the 49 circular sections are addi¬ 
tionally connected via the potential 0 in order to obtain stable 
circles. The two pole particles are only bonded to their 60 near¬ 
est neighbors. The bond strength = 10 4 ksTla 1 is chosen for 
all bonds, whereas the preferred length r e of each bond is deter¬ 
mined by the geometry. Here, kg is the Boltzmann constant and 
T the temperature. 

A flagellum is described by the helical wormlike chain 
model, 13 57 with an adaptation suitable for the combination with 
MPC. As shown in Fig.|2][b), a helical flagellum consists of N = 76 
segments with a total of 381 particles. In each segment, six parti¬ 
cles are arranged in an octahedron of edge length a/\fl, forming 
12 bonds along the edges and three along the diagonals. The pre¬ 
ferred bond lengths are r e = a/\fl for edges and r e = a for diag- 



6 flag., random 8 flag., symmetric 


Fig. 1 Model of E. coli. Bacteria differ in the number of flagella and 
arrangement of flagella on the cell body. In a ‘symmetric’ arrangement, 
the first contour particle of each flagellum is uniformly distributed along 
a circle on the body. In a ‘random’ arrangement, the first contour particle 
is randomly located on the body. 



(b) 
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Fig. 2 (a) Model of the spherocylindrical cell body of diameter d = 0.9 
pm and length 4 = 2.5 pm. It is composed of 51 circular sections of 
particles, which are connected by the bond potential of Eq. (TJ. (b) The 
flagellum, a three-turn left-handed helix of radius R = 0.2 pm, pitch 
A = 2.2 pm and contour length L c = 7.6 pm (corresponding to the parallel 
length Ly = 6.6 pm), consists of 76 consecutive segments. 


onals. The octahedron construction allows for a straightforward 
description of the intrinsic twist of the flagellum and a coupling 
of the twist to the forces exerted by the MPC fluid. 

The bonds b 3 = r i n+l — r i n with n = 1, ...,2V specify the contour 
of the flagellum, and, together with b* = r in+ \ — r in+ 3 and b 7 = 
r in+ 2 — Y i n + 4 , define orthonormal triads {e 4 , e 7 , e 3 }, where ejf = 
bn/|b{?| with a e {1,2,3}. Here, the r i n denote the positions of 
the backbone particles, and the r in+k (k = 1, 2, 3, 4) the positions 
of the particles in the plane with the normal e 3 . 

To characterize the local elastic deformation of a flagellum, 
the triad {e 4 , e 7 , e 3 } is transported to {e* +1 , e 7 +1 , e 3 +1 } along the 
chain. This process is performed in two steps: (i) the rotation 
of {e 4 , e 7 , e 3 } around e 3 by a twist angle %, and (ii) the rotation 
of the twisted triad {e 4 , e 7 , e 3 } by a bending angle around the 









normal n n = (e 3 x e 3 +1 )/|e 3 x e 3 +1 | to the plane defined by the 
contour bonds b 3 and b 3 +1 . The elastic deformation energy is 
then 

u* = l t K ?i l'( n “ - n e ) 2 . (2) 

z a= 1 n= 1 

where is the bending strength, the twist strength, 

and Cl n = fl^e 1 + £!„€% + £2 3 e 3 = strain vector. We 

choose Kl x = K ^ = K^ = 5 • 10 4 &#r, corresponding to a bending 
stiffness of 2 • 10 -23 N m 2 for flagellar filaments within the experi¬ 
mental range of about 10 -24 — 10 -21 N m 2 . The parameters Q® in 
Eq. ([ 2 ]) define the equilibrium geometry of the model flagellum 
and are chosen to recover the shape of an E. coli flagellum in the 
normal state, i.e., a three-turn left-handed helix of radius 0.2 pm 
and pitch 2.2 pm.QSlsee the Appendix | a| for details. 

A flagellum is attached to the cell body by choosing a body par¬ 
ticle as its first contour particle ( 4 , see Fig.[2][b) for the notation). 
The rotation of the flagellum is induced by a motor torque T de¬ 
composed into a force couple F and —F acting on particles i\ +2 
and i\ + 4 (T = b 2 x F with F antiparallel to b{), or equivalently 
i\ +1 and q + 3 (T = b}xF with F parallel to b 2 ). Hence, there is 
no net force on the bacterium. We consider |T| < lOOO/^T ~ 4100 
pN nm, smaller than the stall torque of approximately 4500 pN 
nm of the flagellar motor.ESAn opposite torque —T is applied to 
the body to ensure that the bacterium is torque-free. We do not 
explicitly model the hook that connects a flagellum and the body, 
in contrast to the model of Ref. 20, but incorporate the physical 
features of the hook: (i) transmitting the motor torque to the flag¬ 
ellum for rotation, and (ii) provide the flagellum the freedom to 
adopt any orientation relative to the body. To prevent a flagel¬ 
lum from crossing the cell body or another flagellum, we use the 
repulsive Lennard-Jones potential 

4e[(cx/r) 12 -(cx/r) 6 ]+£, r<^2 a ^ 

0 , otherwise 

with £ = kfiT to capture excluded-volume interactions. For 
flagellum-body repulsion, r is the distance between a flagellar 
contour particle and a body particle, and a = 0.5 a is equal to the 
section spacing of the body. For flagellum-flagellum repulsion, r 
is the closest distance between contour bonds of two flagella, and 
a = 0.25 a is set by the filament diameter around 25 nmJSl 

The dynamics of the bacterium is determined by the forces re¬ 
sulting from the potentials in Eqs. 0-0 and the forces for gen¬ 
erating the torques T and —T, and by the momentum exchange 
with MPC fluid via the collision rule described by Eq. 0 below. 

2.2 Fluid model: Multiparticle collision dynamics 

The MPC fluid is modeled by a collection of point-like particles of 
mass m. Their dynamics proceeds by alternating streaming and 
collision steps. In the streaming step, the fluid particles move 
ballistically and the position r* of particle i with its velocity v* is 
updated according to 

r i(t + At) = 17 (f) + v,-(f) At, (4) 


where At is the time interval between collisions. In the collision 
step, all particles are sorted into cubic cells of length a and the 
velocity v* of particle i in cell c is renewed via the collision ruleEil 

v?ew = Vc + v ran _ £ jec m^ jec mj 

+ - r c ) x (v,- - v) an )] x (r ; - r c ), (5) 

where v c and r c are the center-of-mass velocity and position of 
the particles in c, mj the mass of particle j in c, v^ an a random 
velocity sampled from the Maxwell-Boltzmann distribution, and 
I the moment-of-inertia tensor of all particles in c. The collision 
rule (i) conserves both linear and angular momentum in each 
cubic cell, (ii) includes thermal fluctuations of the fluid, and (iii) 
maintains a constant temperature. To satisfy Galilean invariance, 
a random shift of the collision-cell grid is performed before each 
collision step ED 

2.3 Parameters 

The length a of a collision cell, the mass m of a MPC particle, and 
the thermal energy IcbT define the length, mass, and energy units 
in our simulations. Other units are easily derived, e.g., density 
Po = a~ 3 , velocity vo = \JkgT /m, time to = a^/m/ksT, and shear 
viscosity qo = V m kgT/a 1 . 

The simulations are performed in cubic boxes of length up to 
L = 250 a with periodic boundary conditions. The largest system 
contains more than 1.5 x 10 8 fluid particles. We choose the colli¬ 
sion time step At = 0.05^o and the fluid density p = 10po, lead¬ 
ing to the fluid viscosity q = 7.15 t7q an d the Schmidt number 
Sc = 20, for which the fluid exhibits liquid-like dynamics . 61 New¬ 
ton’s equations of motion for the bacterium model are integrated 
with the time step 8t = At/25 using the velocity-Verlet algorithm. 
The Reynolds number Re = dpv/rf is within the range 0.01-0.2 for 
the simulated bacteria with body diameter d = 9a and swimming 
speed v = 0.00075 — 0.0125vo. 

Matching the geometry of our model bacteria to E. coli with 
a body of length of 2.5 pm and diameter 0.9 pm3 leads to the 
physical length a ~ 0.1 pm of a collision cell. The thermal energy 
is ksT ~ 4.1 pN nm at T = 300 K. A comparison with the viscosity 
of water of approximately 9 x 10 -4 N s/m 2 with rj = 7.15 qo of the 
MPC fluid gives the physical time scale to ~ 30 ps. 

3 Results and discussion 

3.1 Hydrodynamic friction coefficients of the bacterium 

3.1.1 Flagellum 

The flagellar hydrodynamic coefficients relate the force & and 
torque exerted by the fluid on the flagellar bundle to its trans¬ 
lational velocity v and rotation rate co (angular velocity — co by 
right-hand rule) via 

«^ = -^v + of(-to), (6) 

& = +7fv-tf(-co), (7) 

where yf and yf are the translational and rotational friction co¬ 
efficients about the helical axis, and y^ is the rotation-translation 
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Fig. 3 Results for a single flagellum from simulations in cubic boxes of length L = 12 pm: (a) translational speed v versus axial force and (b) 
confining torque % Q nf versus v for a translating nonrotating flagellum; (c) rotation rate co versus axial torque <% ot and (d) confining force ^ coa f versus co 
for a rotating nontranslating flagellum. The red solid arrows indicate the direction of translation or rotation caused by & w \\ or # ot , whereas the red 
dashed arrows imply the direction of rotation or translation which would have arisen from the rotation-translation coupling, but are prevented by ^ conf 
or J^conf from the harmonic potentials on the flagellum. 


coupling coefficient. These coefficients, known as elements of the 
propulsion matrix,^ are essential quantities of bacterial swim¬ 
ming. We determine these coefficients by simulating a translat¬ 
ing non-rotating (co = 0, 0) flagellum and a rotating non¬ 

translating (v = 0, co Y 0) flagellum in the absence of the cell 
body. In the first case, the flagellum is pulled by a force & w \\ 
along its helical axis, say v-axis, and y-, z-positions of each of 
its contour particles are trapped by harmonic potentials. This 
yields a confining torque ^ 0 nf that prevents the rotation by coun¬ 
teracting the torque in Eq. 0- In the second case, a flag¬ 
ellum is rotated by a motor torque, and harmonic potentials 
are applied on its contour particles that prevent its translation, 
but allow for the rotation about the helical axis, causing a net 
axial torque ^ ot and a confining force ^ C onf opposite to the 
force in Eq. (J6j) . The two measurements lead to ^ = ^pull/ v > 
if = -^onf/v, Yf = -^rot/co, and Y = ^confM In the latter 
measurement, two no-slip wallsESSI orthogonal to the flagellum 
are included to ensure that the flagellar speed v vanishes with 
respect to the background fluid at rest. Without such no-slip 
walls, the translational speed of the flagellum is still v relative 
to the fluid velocity at the boundaries since the rotating flagellum 
is constantly pumping fluid backwards. From the results shown 
in Fig. [ 3 ] we obtain y^/(m/to) = 1028 ±2, yf/y/rnksT = -382 ±7, 
Hf/iksTto) = 11051 ±89, and yf/y/rnksT = -400 ± 1. The two es¬ 
timates of yf agree within 5%, confirming the consistency of our 
measurements. 


3.1.2 Cell body 

In a similar way, we measure the translational and rotational fric¬ 
tion coefficients and of the cell body by considering the body 
under an axial force ^ w \\ or torque ^ ot . We obtain y£ = ^ w \\/v = 
(1118±5)m/fo and = ^ G t /co = (52435 ±236)&#77 0 from simu¬ 
lations with box size L= 12 pm. These values are consistent with 
the theoretical values y^/(m/to) = 1023 and YJ(k B Tto) = 46618 
for a solid spherocylinder of the same size J 65 ! 6 ^ The quantitative 
agreement confirms that the applied particle-based mesoscale hy¬ 
drodynamics model can well capture the no-slip boundary condi¬ 
tions on the bacterium surface, as already demonstrated for hard 
sphere colloids in Ref. 67 

Table 1 Comparison of hydrodynamic properties for the model 
bacterium with experimental results of E. coli^f and / are the 
translational and rotational friction coefficients, and f is the 
rotation-translation coupling coefficient. The subscripts ‘b’ and T refer to 
body and flagellum, respectively. All the coefficients are rescaled by the 
absolute values of 7 J, which are 1281 and 1118 m/to from simulations in 
cubic boxes of lengths L = 8 and 12 pm, and yfc = 1.4 x 10~ 2 pN s/pm 
from experiments. 
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3.2 Comparison with experimental results and resistive- 
force theory 

Table [l] summarizes the hydrodynamic coefficients of the bac¬ 
terium body and a flagellum. The values measured for the sys¬ 
tem sizes L = 8 and 12 pm differ by at most 10%, implying that 
the periodic boundaries in the simulations have a rather small ef¬ 
fect. The experimental values in Table [l] are from measurements 
of swimming E. coZi® with an average body length of 3.0 ±0.8 
pm, close to 2.5 pm in our model. The flagellar properties deter¬ 
mined from experiments are for a bundle with, on average, 3 — 4 
flagella, rather than for a single flagellum. In addition, the cell 
body is found to wobble around the swimming axis. Given these 
caveats, the hydrodynamic properties of our model flagellum and 
cell body are in good agreement with experimental results for E. 
coli. 

Table 2 Comparison of the friction coefficients of a flagellum with 
estimates by resistive-force theory (RTF)® and experimental results of 
E. coliJzHjhe theoretical expressions derived by resistive-force theory 
are presented in App.[C] For the current model, only results for the box 
size L = 12 pm are presented. 
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Table [2] compares the friction coefficients with results obtained 
by resistive-force theory (RTF) theoryi^ ^ 25 l 68 l we use the theo¬ 
retical expressions derived in Ref.;68j the terms for the respective 
friction coefficients are summarized in App. [c] In order to avoid 
ambiguities by different expressions for the factors K n and K t of 
Eqs. ( [2l] ) and ( [22] ), as presented in Refs. I2H25I68I we discuss 
ratios of the various friction coefficients only. For the theoretical 
results, we use the pitch angle £ = 0.519, the pitch length A = 2.2 
pm, and the hydrodynamic radius r = 0.1 pm of the flagellum 
string, which yields % = 0.68 (cf. App. |cj). The ratio of/yf is in 
close agreement for the various approaches. However, ratios for 
the simulation model including are about twice larger than the 
prediction of RTF or even three times compared to those deter¬ 
mined from experimental results. As far as the comparison with 
RTF is concerned, the agreement is reasonable, considering the 
fact that the simulation model uses a discrete representation of 
the helix, and the ambiguities in the various parameters of the 
analytical approach. Here, a more thorough comparison between 
theory and simulation, and between the various theoretical ap¬ 
proaches would be desirable. In experiments, the helical bundle 
seems to yield a much smaller rotational fiction coefficient. To 
which extend this friction coefficient is affected by the nature of 
the bundle and the above mentioned wobbling dynamics needs to 
be studied further. 

3.3 Swimming speed and flagellar rotation rate 

In Fig.[4ja), the bacterial swimming speed v is displayed as a func¬ 
tion of the rotation rate co of the flagellar bundle for our model 
bacteria of Fig. [l] which differ in the number of flagella and the 
arrangement of flagella on the cell body, ranging from one to eight 


flagella with symmetric and random arrangements of the anchor¬ 
ing points. The rotation rate of the flagellar bundle is co = 2tz/z, 
where z is the average time for the helical bundle to complete one 
revolution. The data points in Fig. |4fa) are obtained from exten¬ 
sive hydrodynamics simulations, where each bacterium travels a 
distance larger than its own length of about 10 pm. The solid 
line in Fig. [ 4 ] is a least-square fit to all data points through the 
origin, yielding v/coR = 0.07, comparable to the experimental ra¬ 
tio v/coR = 0.14,® where R = 0.2 pm is the helix radius of the E. 
coli flagellum. The scatter of the data points around the linear 
relation may arise from the fact that individual bacteria exhibit 
different wobbling amplitudes of the body. 



Fig. 4 (a) Bacterial swimming speeds v as function of the rotation 
frequency co of the flagellar bundle and (b) rotation frequency co as 
function of the flagellar motor torque |T| for the model bacteria shown in 
Fig.Q The solid line in (a) is a least-square fit to all data points. The 
solid lines in (b) are least-square fits to the data points in the same color. 

The ratio of swimming speed v to flagellar rotation rate co 
can be understood from the hydrodynamic properties of the bac¬ 
terium. The force in Eq. © is balanced by the translational 
friction force — exerted by the fluid on the cell body, which 
gives 

v 

- = ( 8 ) 
co f h + i f 

We obtain v/coR ~ 0.1 and 0.13 using the model and experi- 













mental values in Table [T] respectively. These ratios are close to 
v/coR = 0.07 from the fitted line in Fig. [4|[a) and the experimen¬ 
tal ratio v/coR = 0.14. The agreement emphasizes the importance 
of hydrodynamic coefficients for a quantitative understanding of 
bacterial swimming. 

Figure |4|[b) shows that the rotation rate of a bundle formed by 
four flagella near one pole of the cell body (red points) is larger 
than that of a single flagellum (black points) at the same motor 
torque, e.g., co/In = 340 Hz versus 223 Hz at |T| = 500ksT. This 
result is consistent with the experimental observation ® that a 
flagellar bundle rotates faster than a single flagellum. The en¬ 
hanced rotation as a result of the hydrodynamic coupling be¬ 
tween the flagella and has also been found in simulations of he¬ 
lical bundles,® where one of the ends of the helices is fixed in 
space. We note, however, that such hydrodynamic enhancement 
can be outweighed by the flagellum-body friction, depending on 
the connection of flagellar filaments to the body. As shown in 
Fig.gb), the bundle of four flagella with a ‘symmetric’ arrange¬ 
ment on the body (purple points, Fig. []J d)), rotates slower than 
both a single flagellum (black points, Fig. [lja)) and the bundle 
of four flagella near one pole of the body (red points, Fig. [ljc)) 
at the same motor torque |T| < 400 ksT. In E. coli, each flagel¬ 
lum is attached to the body via a flexible hook of approximately 
50 nm length,® much shorter than the flagellar helix radius of 
0.2 pm . The rotating flagella that form a bundle spanning over 
the body are very likely to collide with the body, contributing to 
the flagellum-body friction. With the single-flagellum rotational 
friction coefficient = 11051&g77o, the torque rotating the bun¬ 
dle of four flagella with co/2n = 196 Hz (co = 0.037 /to, see the 
purple point at |T| = 500 kgT) is rfco ~ 409 hgT, only about one- 
fifth of the total motor torque 2000 k^T, indicating a significant 
flagellum-flagellum and flagellum-body friction. 

3.4 Flow field 

Figure [5] shows the time-averaged flow field created by a model 
bacterium of approximately 8 [im length in a periodic cubic simu¬ 
lation box of length L = 25 p m. Each of the four flagella is turned 
by a motor torque |T| = 1000 kgT, leading to a bundle with ro¬ 
tation frequency co = 0.074/to. The bacterial swimming speed is 
v ~ 0.0125vq and the propulsion force / p = ~ 0.57 pN consis¬ 

tent with the experimental values of about 0.4-0.6 pN. 2 ®^ Har¬ 
monic potentials are applied to y- and z-positions of the center- 
of-mass of the body and flagellar bundle such that the swimming 
axis is parallel to the v-axis. The flow field in the swimming plane, 
shown in Fig.[5ja), is computed by averaging the velocities of the 
fluid particles cylindrically symmetrically around the swimming 
axis. The flow fields in Figs. [5jc)-(g) are obtained as time av¬ 
erage of the fluid-particle velocities in the planes perpendicular 
to the swimming axis at the different locations indicated by the 
white vertical lines in Fig.[5ja). 

The flow pattern not too close to the bacterium approximately 
resembles that of swimming E. coli determined from experiments 
(see Fig. 1A in Ref. 29). Closer to the bacterium, the flow field 
exhibits specific features reflecting the bacterium’s detailed struc¬ 
ture. In particular, the flow field reveals a front-back asymmetry, 


since the cell body and flagellar bundle are physically different 
units. Along the cell axis, we find low fluid velocities in front of 
the cell, behind the cell body as well as in the middle and behind 
the flagellar bundle. The flow velocity is high along the whole 
cell body, at the flagellar bundle somewhat behind the cell body, 
and toward the end of the bundle. The streamlines are closed 
in Fig. [5](a) as a consequence of the applied periodic boundary 
conditions, which implies differences in the far field compared to 
experimental observations. The effect of the boundary conditions 
is confirmed by the theoretical flow field for a finite-distance force 
dipole with the same boundary condition in Fig.[5fb). 

The flow patterns in the planes perpendicular to the swimming 
axis illustrate the interplay between the rotating flagellar bundle 
and counterrotating cell body. As shown in Figs. [5jc), (e), and 
(g), the flow field exhibits two spiral vortices associated with the 
rotation of the cell body and flagellar bundle, respectively. In 
Fig-|5fd), the fluid follows the clockwise rotation of the flagellar 
bundle in the central region, and tends to rotate counterclockwise 
in the outer region, implying a strong influence of the cell body. 
In Fig. [5jf), the fluid rotates counterclockwise as the body and 
there is no significant effect from the flagellar rotation. 

The flow field of the discrete-particle bacterium model of Ref. 
l20l with a triangular-prism body and three flagella, each com¬ 
posed of 15 beads, involves an infinite fluid domain, because hy¬ 
drodynamic interactions are captured by the Rotne-Prager hydro- 
dynamic tensor. Consistent with this study, we find spiral and 
helical flows (Figs. [5fc) -(g)). However, there are also distinct dif¬ 
ferences, namely (i) three vortices are predicted in Fig. 5 of Ref. 
l20l and (ii) fluid flows toward the flagellar bundle from behind 
(Figs. 4 and 7) rather than away from the cell as in Fig. [5] It is 
not a prior evident where these differences come from. A reason 
could be the more detailed modelling of the cylindrical cell body 
and the flagella in our approach. Furthermore, in agreement with 
our studies, the experimentally determined flow field presented in 
Ref. [29] shows no indication of a forward flow at the end of the 
flagellar bundle. 

We analyse next, how well the flow field can be reproduced 
by a simplified description in terms of a force-dipole and rotlet- 
dipole model. Figure |6j[a) displays the fluid velocity along the 
swimming axis both in front of the body and behind the flagellar 
bundle. The theoretical result for the force dipole, with the force 
/ p ~ 0.57 pN and the 4 = 3.84 ^m (see App.jBj), agrees quite well 
with the simulation data for distance r/4 > 1 , he., more than one 
body length away from the cell. The flow for r/4 > 3 decays much 
faster than r~ 2 due to the periodic boundary conditions. In this 
regime, the flow field is consistent with the result for the finite- 
distance force dipole of two Stokeslets with periodic boundary 
conditions (cf. App. [b]). In the range 0.8 < r/4 < 1-8, Vfl u i d de¬ 
creases approximately as r -3 , a dependence previously observed 
in Ref. [20 Already for r/4 > 2, our flow field agrees well with 
that of the finite-distance force dipole. Hence, up to r/4 <3 — 4 
the flow field is hardly affected by the boundary conditions. The 
asymptotic force-dipole dependence Vfl u i d ~ r -2 is only reached 
for distances r/4 > 5. Figure |6][b) displays the flow field in 
planes perpendicular to the swimming axis through the center 
of mass of the body or flagellar bundle, respectively. The Simula- 



Fig. 5 Time-averaged flow field generated by a single swimming bacterium as obtained from simulations: (a) flow field in the swimming plane (b) the 
theoretical flow pattern for a finite-distance force dipole as illustrated in Fig.[6ja) as superposition of two Stokeslets within the same periodic box as for 
our simulations), (c) -(g) Flow fields in planes perpendicular to the swimming plane at positions indicated by the white vertical lines in (a). The 
streamlines indicate the flow direction, and the logarithmic color scheme indicates the magnitude of the flow speed scaled by the bacterial swimming 
velocity. 


tion results for r/£\>> 1 are well described by a superposition of 
a force dipole with force f p ~ 0.57pN and a rotlet dipole with the 
torque N ~ 80one-tenth of the torque rotating the flagellar 
bundle y£&> ~ 800&g7\ The rotlet-dipole contribution dominates 
at r/4 < 2, whereas the force-dipole contribution dominates at 
larger distances. This suggests that the counterrotation of cell 
body and flagella needs to be taken into account for an appro¬ 
priate description of near-field hydrodynamics of bacteria. Our 
result provide the first measurement of the rotlet-dipole strength, 
which is important for a quantitative description of near-surface 
swimming behaviors of bacteriaJSED it would be interesting to 
compare our predictions with future experimentally determined 
rotlet-dipole strengths for flagellated bacteria. 

4 Conclusions 

We have developed a coarse-grained model for an E. coli- type 
bacterium and investigated its swimming behavior by mesoscale 
hydrodynamics simulations. The suggested implementation of 
octahedron-type segments to construct a model flagellum of he¬ 
lical structure with bending and torsion elastic energy is particu¬ 
larly suitable in combination with a particle-based fluid, because 
it provides a simple means to capture the twist by the fluid on 
a flagellum. The hydrodynamic friction coefficients of the model 
bacterium and the relation between its swimming speed and flag¬ 


ellar rotation rate both quantitatively agree with experimental re¬ 
sults for E. coli. The flow field created by the model bacterium 
exhibits a rather complex pattern adjacent to the bacterium, with 
and a two-vortex structure due to the counterrotation of the cell 
body and flagellar bundle. From comparison to theoretical pre¬ 
dictions for force dipole and rotlet dipole, we find the simplified 
dipole model can well describe the flow more than one cell-body 
length away from the bacterium, and the rotation of the body 
and flagella makes a dominant contribution to the near-field flow. 
We provide the first measurement of the rotlet-dipole strength, an 
essential quantity for modeling bacterial motion near surfaces. 

The simulations show that our model is very well suited for 
theoretical studies of swimming bacteria. By adjusting the geo¬ 
metric parameters, one expects to achieve a similar quantitative 
description for other bacteria such as Bacillus subtilis, Saolmonella 
typhimurium, Rhodobacter spheroides, and Rhizobium lupini. 

Naturally, further aspects of swimming bacteria can be investi¬ 
gated by our approach. By modifications of the elastic deforma¬ 
tion energy in Eq. © to account for flagellar polymorphic trans¬ 
formations EE0 and by implementing reversal of the flagellar mo¬ 
tor torque T, the full run-and-tumble motion can be addressed. 
Since boundaries and external flow are easily implemented in the 
MPC method, our simulation approach opens an avenue for de¬ 
tailed studies of confinement effects 51 and non-equilibrium as- 


















































Fig. 6 The fluid velocity v fluid as a function of distance r to the center of 
mass of the body or flagellar bundle (a) along the swimming axis, and 
(b) in the planes perpendicular to the swimming axis through the body or 
flagellar center of mass, respectively, v is the swimming speed and 4 
the body length. In (a), the solid line is the theoretical result for a force 
dipole of strength / p 4, where / p ~ 0.57 pN is the bacterial propulsion 
force and 4 = 3.84 pm the center-of-mass separation between the body 
and flagellar bundle. The dashed line is for the force dipole with the 
same periodic boundary conditions (PBC) as in simulations. In (b), the 
red and green solid lines represent the results for the force dipole as in 
(a) and the rotlet dipole with torque N ~ 80 k B T, one-tenth of the torque 
rotating the flagellar bundle co ~ 800& 5 7\ The cyan solid line is the 
superposition of the force dipole and rotlet dipole. 


pects in bacteria locomotion. 

A Flagellum model 

We estimate the bending and twist strength of our discrete 
flagellum by mapping the elastic energy in Eq. © to a continuous 
form. We identify the directions e 3 of the bond vectors b 3 as local 
tangent vectors on the filament. Then, the discrete measure of 
the local curvature k = |e 3 +1 — e 3 |/<z = 2sin($„/2 )la turns into k « 
D n /a for small bending angles and the local torsion becomes 
t = — |n n+ i — n„| /a = -2sm(\]/ n /2)/a « -\j/ n /a for small torsional 
angles i jf n between the normal vectors n n and n n+ i, where n n = 
(e 3 x e 3 +1 )/|e 3 x e 3 +1 1. At small \\f n ~ (p n since the deformation 
mainly originates from twist. In Eq. the components of the 


strain vector are 


n l ~A 2 3 


sin 


n ' c n+15 


n£ = 


A 

sin A 


c n ‘ c n+15 


&n = <Pn- 


(9) 

(10) 

(ID 


Using e^ +1 = sin A cos <j> n e h + s i n A sin A e n + cos A^n> where A is 
the azimuthal angle of e^ +1 in the plane defined by e, 1 , and ej n 
expansion of Eq. © yields 

1 n-i r i 

i/el = T {4 4 2 + 4 2 -2d„4cos(0„-<4 +4((p„-<p e ) 2 } 

n= 1 


* [«d(0» - A) 2 + 4(<P« - <Pe) 2 ] 

n= 1 L 

[ K ll a ( K ~ K e) 2 + K el a (^~ ^e) 2 ] ds, (12) 

with the of the linearized expressions of Eqs. 0-GD and 
for small deviations of (j) n from its equilibrium value <j> e . The in¬ 
tegral over the flagellar contour describes the continuous elastic 
energy. 57 

The experimental values for the bending stiffness of bacte¬ 
rial flagella range from about 10 -24 to 10 -21 N m 2 .USEE! For 
S. typhimurium, K^a is estimated as 2 — 4 x 10 -24 N m 2 from 
quasi-elastic light scattering in Ref. (m consistent with the value 
3.5 x 10 -24 N m 2 measured from pulling experiments in Ref. 
[161 In the Ref. [72] a Young’s modulus of about 10“ 11 N m -2 
has been determined for S. typhimurium and E. coli from elec¬ 
tron micrographs, which yields the approximate bending stiff¬ 
ness 10 -21 N m 2 given a filament radius of 10 nm. The twist 
or torsional stiffness has been generally assumed to be of the 
same order of magnitude as the bending stiffness. Our choice 
= K= 5 x 10 4 IcbT, with the physical length scale am 0.1 
pm, corresponds to a bending stiffness of K^a ~2x 10 -23 N m 2 
well within the experimental range. 

In order to determine the equilibrium parameters in Eq. 
we consider a flagellum in the normal state—a three-turn 
left-handed helix of radius 0.2 ^m and pitch 2.2 ^m. The 
bending angle between subsequent discrete bond vectors is then 
$ e = arccos (e 3 • e 3 +1 ) = 0.125 for the bond length a = 0.1pm. We 
choose e* in each segment such that it coincides with n n after the 
twist with angle cp G , i.e., e* = n n . e* is thus invariant in bend¬ 
ing from e 3 to e 3 +1 around n n by the angle hence, e* +1 = 
e* = e 3 x e n-uJ sint^e- Then, <p e = -arccos(e* • e* +1 ) = -0.217. 
From Eqs. Q9])— (flT]), we obtain fig = 0.122, fl 2 = —0.027, and 
fl 3 = —0.217, which are the values used in our simulations. 

B Force-dipole flow field for periodic 
boundary conditions 

We consider the flow generated by a force dipole as illustrated by 
the cartoon in Fig.^a). The two point forces fi = (—/ p ,0,0) and 
f 2 = (/ P ,0,0) ^ located at ri = (-4/2,0,0) and r 2 = (4/2,0,0), 














respectively. The fluid velocity at position r is given by 
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u(r) = G(r-ri)-fi+G(r-r 2 )-f 2 , (13) 


with the Oseen tensor 


G(r) = 


1 / r0r\ 

8ttt 7 |r| V + |r| 2 / 


(14) 


To calculate the flow of the force dipole in a cubic box of length L 
with periodic boundaries, one can start with the Oseen tensor in 
k-space 


G(k) 


1 / k(g)k\ 

T7]kp V |k|^7 * 


(15) 


The k-space fluid velocity is then 


u(k) = G(k)f!(k) + G(k)f 2 (k), (16) 


with the Fourier-transformed point force f) (k) = /f;<5(x — 
Ti)e~ lk ' x d 3 x = fie~ ± ri , i = 1,2. Inverse Fourier transformation for 
the periodic system yields the flow field 

u'(r) = £ u(ky k ' r (17) 

k^O 


with the three components of k assuming multiples of 27 i/L. 


C Resistive-force theory for helix 


Expressions for the friction coefficients of Eqs. ([§]) and ([7]) for 
a helix have been presented in Refs. 21 22 25 68 The original 
expressions of Ref.f68lare 


(18) 


(19) 


(20) 


with % = K t /K n , £ = tan where £ is the pitch angle, and 


4 - 717 ] 

ln(cA/r) + 1/2’ 
270 ] 

ln(cA/r) 


( 21 ) 

( 22 ) 


Here, A is the pitch, r the (hydrodynamic) radius of the flagella 
filament, and c = 0.18 a constant.®] For % = 1/2, which corre¬ 
sponds to the ratio of the viscous coefficients of an infinitely long 
rod, and c = 2, the Eqs. ( [18] ) - ( [22] ) reduce to the definitions in 
RefJ25 The viscous coefficients of Ref. 21 are somewhat different, 
specifically the expression for K t . 
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